p-wave phase shift and scattering length of 6 Li 



S. GautarrQ and D. Ang 

Physical Research Laboratory 
Navarangpura, Ahmedabad - 380 009 

We have calculated the p-wave phase shifts and scattering length of 6 Li. For this we solve the p 
partial wave Schrodinger equation and analyze the validity of adopting the semiclassical solution to 
evaluate the constant factors in the solution. Unlike in the s wave case, the semiclassical solution 
does not provide unique value of the constants. We suggest an approximate analytic solution, which 
provides reliable results in special cases. Further more, we also use the variable phase method to 
evaluate the phase shifts. The p-wave scattering lengths of 132 Cs and 134 Cs are calculated to validate 
the schemes followed. Based on our calculations, the value of the p wave scattering length of 6 Li is 
-45a D . 

PACS numbers: 34.50.-s,34.10.+x 



I. INTRODUCTION 

For bosonic isotopes of atoms, one parameter which 
describes the low energy scattering properties is ao, the 
s-wavc scattering length. It arises from the s-wave phase 
shift rj . The corresponding interaction potential is a 
crucial parameter in Bose Einstein condensates of di- 
lute ultracold atomic gases. The fcrmionic counterpart is 
the interaction potential arising from the p-wave scatter- 
ing. Calculations of which is important as experiments 
on fermionic isotopes have made impressive strides since 
the first experimental observation of degenerate fermions 
1]. Superfluidity in two species fermionic mixture of 6 Li, 
which was first predicted theoretically 0], has been ob- 
served without any ambiguity 0, 0, Of - Following which, 
intense experimental and theoretical investigations con- 
tinues on the phase diagram of the spin polarized two 
component 6 Li mixture. Among the most recent devel- 
opments are the theoretical investigation of the phase 
diagram at finite temperature [f| and the experimen- 
tal investigation of the same at unitarity Q. Further, 
the recent achievement of cooling the 6 Li- 40 K fermionic 
mixture to degeneracy Q takes us closer to observing 
exotic phases predicted for spin polarized hetcronuclcar 
fermionic mixtures Q . Some of the predicted phases are 
fragile and crucially dependent on the difference of the 
chemical potentials. For such cases, the change in chemi- 
cal potential induced by the p-wave scattering is likely to 
be an important parameter. The pscudopotcntials aris- 
ing from the higher partial waves, discussed in recent 
works [l(| [H| j can be used to incorporate the effects of 
p-wave scattering in fcrmionic isotopes. 

In this paper we describe the calculation of the p-wave 
scattering length of 6 Li. As suggested in an earlier work 
on s-wave scattering [l2j], the WKB method is used to 
determine the constant parameters of the partial wave 
solutions of the Schrodinger equation. However, unlike 



in the s-wave scattering, the parameters in the p-wave 
calculations has radial dependence. This is an outcome 
of including the centrifugal potential in the effective in- 
teratomic potential. To circumvent this, an analytic ex- 
pression, which provide an estimate of the p-wave phase 
shift is suggested. This method is valid when the disper- 
sion constant Cq is large. We also calculate the phase 
shift using the variable phase method. To test and val- 
idate the numerical schemes adopted, the s-wave phase 
shift of 133 Cs and p-wave scattering length of 132 Cs and 
134 Cs are calculated. 

For completeness in Section. I of the paper, we provide 
an outline of solving the p partial wave Schrodinger equa- 
tion in different radial regions. The Section. II discusses 
the calculation of the scattering length using the WKB 
method to determine the constants in the partial wave 
solutions. Then a brief description of the variable phase 
method is provided in Section. Ill, this is followed with 
results and conclusions. All the equations and results in 
this paper are in atomic units, in which Ti = m e = e = 1. 



II. p-WAVE PHASE SHIFT 

The radial part of the Schrodinger wave equation for 
collisions between two atoms with V(R) as the inter- 
atomic potential is 



dR 2 



k 2 - U(R) - 



1(1 + 1) 
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X = 0, 
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Here U(R) = 2mV(R), k is the relative momentum of the 
two atoms, I is the angular momentum quantum number 
and is the radial wave-function. This equation can 
be used to calculate scattering phase shifts (iji) which 
each component angular momentum suffers due to the 
interaction with the scattering center. At low energies 
we can calculate approximate solution of Eq. ([I]) in three 
difference radial ranges. These are the kR <g; 1, kR < 1 
and kR 1 regions. 
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A. kR <C 1 region 



C. kR ~ 1 region 



In this region, we can neglect the k 2 term 
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R 2 



x = o, 



(2) 



and hence, the solution is independent of k. At large 
distances from the scattering center, van der Waal's po- 
tential U(R) = -2ma/(R 6 ) = ~^ 2 /R 6 is the dominant 
inter-atomic interaction. Here, a is the van der Waal's 
coefficient. For low energy collisions k — > 0, the inter- 
atomic potential approaches the asymptotic form within 
the kR <C 1 region. Then, 
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Substituting X (R) = £(R)V~R and R = VtT 2 ^, wc ob- 
tain the Bessel differential equation of order (21 + l)/4 
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The general solution of Eq.([3]) is 
X (R) = VR 
Considering the p(l = 1) partial wave 



= 0. 



(^2 ) + BJ -(2l+l)/i (2^) 



(4) 



- 3/4 (2M 

In the limit i? — * 00 or x — > 0, where fci? < 1, the leading 
terms in the series expansion of J^/^x) and J-3/a( x ) are 
dominant and the remaining terms are negligible 



x(R) = B 



( 7 /2)- 3 / 4 i? 2 



A- 



(7/2) 
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(5) 



2- 3 / 4 r(i/4) 2 3 / 4 r(7/4)i?' 
where T(. . .) are the gamma functions. 

B. kR < 1 region 

In this region, we can neglect the interatomic inter- 
action potential U(R), if it approaches zero faster than 
1/R 2 . In addition, as in the kR 1 region, the k 2 term 
can be neglected 



d 2 X 1(1 + 1) 



x = o. 



(6) 



dR 2 R 2 
The most general solution of this equation is 

X(R) = aR l+1 + bR~ l . (7) 

For 1=1, comparing the solutions in Eq. §5§ and Eq. ([7]) 
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(7/2) 



-3/4 
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and b = A 
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The comparison is possible as x m the R — > 00 of Eq.([3|) 
has the same form as the solution in the fci? < 1 region. 



Increasing R further, we enter the kR w 1 region. In 
this region, we can neglect U(R) in comparison to the 
other terms in Eq. ([T]), then 
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This is the Schrodingcr equation in the asymptotic region 
and the solution is a plane wave of wave number k. How- 
ever, the interatomic potential introduces a phase shift 
to the plane wave 

X (R) = CR 1 ' 2 [J l+1/2 (kR) + (-l) l tsa(rn)J-i-i/2(kR)] . 

Expanding the Bcsscl's functions J l+1 / 2 (kR) and 
J-l-l/iQvR), the solution when kR < 1 is 



x(R) 



t) 



1/2 
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Equating the solutions in Eq. (J7]) and (H 
62 2/ (r(/ + l)) 2 
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Since b/a is independent of k, tan 77; varies as k 2l+1 . It 
should also be mentioned that, though the entire radial 
range is divided into three domains, the solutions are in 
the asymptotic region. This defines the phase shift r/i in 
terms of the b/a and it is strictly applicable for k — > 0. 
For the p-partial wave 



tan rji = — 



or(3)r(4) 



(12) 



In the above expression, to calculate 771, the ratio b/a 
should be determined. From the Eq.©, the phase shift 
can also be defined in terms of A/B as 



tan rji 



1 A 7 3 / 2 r(i/4) 

24 B T(7/4) 



(13) 



Here the actual values of T(3) and T(4) are used. Now 
onwards this definition of rji is used. As the phase shift rji 
has functional dependence on tan -1 , which has a domain 
of [—00, 00], its accuracy is sensitive to the value of A/B. 



III. SCATTERING LENGTH CALCULATION 



As mentioned earlier, the expression of 771 in Eq. (|13[) 
is in the asymptotic domain of the interatomic poten- 
tial and require evaluation of A/B. The solution of the 
Schrodingcr equation Eq.([T]) in the asymptotic region, 
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given by Eq.Q and parametrized in terms of A and B, is 
the outer solution. To determine A/B, the Eq.flT]) should 
be solved within the inner part of the interatomic poten- 
tial: inner wall and well. This can be evaluated using 
the WKB method [T^]. Matching the solutions, outer 
and inner, at a point we can evaluate A/B. The match- 
ing point should be in the region where the asymptotic 
form of the interatomic potential begins to dominate but 
the WKB solution is still valid. 



A. WKB solution 

For partial waves other than s, the centrifugal potential 
is nonzero. Combining the the interatomic potential and 
centrifugal term, the effective interatomic potential 

U eS (R) = U(R)+l(l + l)/R 2 . (14) 

Then the local momentum of the scattered atom p(R) = 
y/—U c ff(R), and the de Broglic wavelength X(R) = 
2 71 "/ yJ~U c ff(R). The WKB approximation is applicable 
when 



(A 



dR 



< 2tt or 



m\F\ 
P 3 



< 1. 



(15) 



In above relation F = —dU e s/(2mdR) is the net force 
acting on the atom. At large distances, when the in- 
teratomic potential approaches the asymptotic form, the 
above inequality for I = 1 is 

3/2 

(16) 

This inequality is in general valid up to large radial dis- 
tances. As an example FigfT] shows that for Cs the in- 
equality is satisfied up to w 30. Then the WKB solution 
at R which satisfies the inequality and larger then the 
classical turning point R Q is 



[3 7 2 


- 2 


<C 


\l 2 1 
J 2 


i? 4 




i? 4 



x(R) = —p co 

VP \JRo 
and the logarithmic derivative is 
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An analytic expression of ((R) for I = 1 is obtained if 
only the van der Waal's potential is considered 
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FIG. 1: Terms in the relation [3 7 2 /i? 4 - 2] < [7 2 /^ 4 - 2 ] 3/2 
which defines the validity of the WKB solution are plotted 
for Cesium. The red and black curves correspond to [y 2 /R 4 — 
2] 3/ ' 2 and [3 7 2 /f? 4 — 2] respectively. Inset shows the ratio of 
the two terms. 



This approximation neglects the well and inner wall part 
of the interatomic potential. It is appropriate for model 
potentials comprising of the van der Waal's and hard core 
potential. The above expression is reduced to that of s 
case when 1 = 0, which is used implicitly in ref. [l2j|. 

From the inequality (|16[) . it is evident that when 7 
is sufficiently large, WKB solution is valid up to radial 
distances where the asymptotic form of the potential be- 
gins to dominate. Choose a point R* in this region, such 
that the 2 in the inequality can be neglected. Then the 
inequality is modified to 



R* < (7/3) 1/2 . 



(20) 



At this point, as the interatomic potential has the 
asymptotic form, Eq. ([4| is a solution of the Schrodingcr 
equation Eq.|T|). This is discussed in Section lll At 



B. Matching x'(R)/x(R) 

From the definition, at R* and around it x = 
7/(2i?* 2 ) » 1. Where x is the variable first used in 
Section IITaI We can then use the asymptotic expression 
expression of the Bessel's function 



lim J n (x) 



cos 



x — \n + 



and the solution in Eq. (j4|) is 
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This is the inner limit of the asymptotic solution. 
Whereas the solution in Eq.®, in the x — > 0, is the outer 
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limit of the asymptotic solution. To determine A/B we 
match the logarithmic derivatives of the analytical solu- 
tion in Eq. j2T]) and WKB solution in Eq. l(T7|) . The ratio 
is 



.4 
B 



\/2ma sin fa - R* 2 [R*({R*) - 3/2] cos0 2 
R* 2 [R*C(R*) - 3/2] cos 0i - V2ma sin fa 



(22) 



here fa = 7 /(2i?* 2 - 5tt/8), 02 = 7/2i?* 2 + tt/8. It is 
evident from the above expression that depends on 
the choice of the matching point R* . This implies that 
r\\ and hence the scattering length, to be defined later, 
depend on R*. This radial dependence arises from the 
centrifugal part in the expression of the effective potential 
U c ff(R). Neglecting the terms arising from the centrifugal 
potential in the expression of C{R)> logarithmic derivative 
of the WKB solution is 



3 

2R' 



(23) 



For the obvious reason, large value of R/j, this expres- 
sion is not valid at very large values of R. Assuming that 
the maximum radial distance at which the above expres- 
sion is valid, the interatomic potential has the asymptotic 
form. If A/B is calculated using Eq. (J23J) as the WKB 
solution, it is then independent of the choice of matching 
point 



.4 
B 



sin — 



7^ 

8 t 



(24) 



where </> — j/(2R 2 ). This is however a crude estimate as 
the outer solution includes the effect of centrifugal term 
but the inner solution, the WKB solution, does not. Still, 
it is a useful estimate, as it is an analytic expression and 
not difficult to evaluate. 



IV. VARIABLE PHASE METHOD 

Another approach to calculate a\ is to evaluate r\\ from 
the phase function equation, which is referred to as the 
variable phase method The method is applicable 

when the interatomic potential satisfies two conditions. 
First condition is, the interatomic potential should be 
less singular at origin than the centrifugal part, that is 



lim \R 2 U(R)} 



0. 



(27) 



This implies that, near the origin U(R) — > UqR™ 1 as R — > 
0, with m > —2, where Uq is a constant. The interatomic 
potential satisfies this condition since at smaller radial 
distances, it approaches the inner wall and is repulsive. 



150 - 



100 - 




100 
R 



FIG. 2: The phase function r)i(R) for 134 Cs at k = 0.05 in 
atomic units. The radial distance is in log scale and the plot 
shows r]i(R) up to R — 500ao, where the the phase function 
converges to an accuracy of 10 -6 . 



C. p-wave scattering length 

For any arbitrary potential which vary as R~ n , the 
scattering length (ai) and phase shift (rji) are related as 

M 



lim k 2l+1 cot f)i 



-1 



fc->0 



,2Z+1 ' 



(25) 



This relation is valid provided I < (n — 3)/2 and is ap- 
plicable for p partial wave when n = 6 or higher. From 
Eq. (Tj"3"l) and (|25[) scattering length of the p partial wave 



1 A 7 3 / 2 r(l/4) 
24B r(7/4) 



(26) 



Using the value of A/B in Eq. (|2"2"l) , the p-wave scattering 
length a\ can be evaluated. A rough estimate of a\ is to 
use the analytical expression in Eq. ([24| . The quantity 
a\ is referred to as the pwave scattering volume. It is an 
important parameter in the three-body recombination of 
identical, spin polarized fermions in low tempcrturcs fl6j . 



The second condition is, the interatomic potential 
should decrease faster than the Coulomb potential or in- 
verse of R. In the variable phase method, the phase shift 
•qi is a solution of the nonlinear differential equation (l3T | 

4(R) = -k- 1 U(R)[coaru(R)ji(kR) - sm m (R) ni {kR)} 2 , 

(28) 

where ji(kR) and ni(kR) are the spherical Bessel and 
Neumann functions respectively. The above equation can 
be used to evaluate the phase function rji{R), which de- 
fines the total phase shift upto R. Then the phase shift 
is the limiting value 



m = lim m{R)- 

R— »OG 



(29) 



To calculate rji, the phase function rji(R) is evaluated 
upto a cut off point at which r]i(R) saturates. The phys- 
ical implication of the cut off point is, it is the radial 
distance beyond which the inter atomic potential can be 
considered zero. To illustrate the radial dependence and 
saturation FigH shows Tfr(R) for k = 0.05 of 134 Cs. The 
figure shows that, for the chosen parameters, saturation 
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occurs at r* 500. The phase shift 77;, as it is evident from 
Eq. I|28p , depends on the relative momentum of the col- 
liding atoms k. To remove the mod(7r) ambiguity, the 
phase shift is normalized such that it approaches zero as 
k — ► 00, that is 

lim m = 0. (30) 

k — >oo 

In addition to above condition 77; is considered to be reg- 
ular function of k. The Eq. is numerically integrated 
from the classical turning point, where the inner wall 
starts, to the cut off point. The solution is the phase 
shift as defined in Eq. ([29]) . 



V. RESULTS 

To evaluate the phase shift and the scattering length, 
the interatomic potential should be known accurately. In 
this paper, we present the results of our calculations for 
cesium and lithium atoms. For which the interatomic po- 
tentials are known accurately and hyperfine interactions 
are less significant in the long range part of the poten- 
tials. The cesium interatomic potential is the one used 
in the work of Gribakin and Flambaum [l2| . which is 
based on ref [13] ■ The lithium interatomic potential is 
based on the work of Zemke and Stwalley . The same 
authors and their collaborators have also calculated the 
interatomic potential of other alkali metals sodium [l9| 
and potassium [2(| • However, these are limited to radial 
distances where hyperfine interactions are not important. 

Once the interatomic potential is known, the loga- 
rithmic derivative of the WKB solution C(-R)i given in 
Eq. (fT8]) is evaluated numerically. To check the accuracy 
and validate the numerical schemes adopted, the numeri- 
cally calculated C(i?) with only the van der Waal's poten- 
tial is compared with the analytic expression in Eq. lflU)) . 
The calculated C(-R) is then used in Eq. (|22]) to evaluate 
A/B. As discussed earlier, the calculated A/B has radial 
dependence and it is difficult to choose an appropriate 
matching point R* . An estimate of R* can however be 
obtained from the results of the variable phase calcula- 
tions. Such a calculation provides a consistency check 
on the use of the WKB method. Despite the radial de- 
pendence, for the p-wave calculations the importance of 
using the WKB solution lies in Eq. ([24|) . an analytic ex- 
pression of A/B. From which, it is possible to calculate a 
rough estimate of the scattering length, which is in rea- 
sonable agreement with the numerical result for 134 Cs. 
For 132 Cs and 6 Li, the two results are very different. 

The phase shift 771 is also calculated from the phase 
function equation, for which the nonlinear differential 
equation Eq. (|28l) is numerically solved . In the present 
work, the Runge-Kutta-Fehlberg (RKF) method is used. 
With this method, it is possible to calculate phase shifts 
for small values of k, which is otherwise not possible 
with methods like fourth order Rungc-Kutta. Calcula- 
tions with the later method has large errors for small 



values of k. For example, in the Cesium calculations 
with RK4 method, the calculated phase shifts are reli- 
able up to k — 6.0 x 10 -3 and k = 1.1 X 10~ 3 for p and 
s partial waves respectively. In comparison, with RKF 
method one can calculate phase shifts for still lower val- 
ues of k. In the numerical calculations, to integrate the 
phase function differential equation, the classical turning 
point is chosen as the starting point. Then, zero is an 
appropriate initial value of r)i(R). Using RKF we can 
calculate 771 (R) for different values of k. The 771 (i?) is 
calculated till a point where it does not increase with 
further increase in R. This point is chosen as the cut 
off point where potential can be considered as zero and 
the corresponding value of the phase function is the re- 
quired phase shift. To calculate the scattering length a\, 
the phase shift 771 is evaluated for a range of k close to 
zero. This is essential as a\ depends on the nature of the 
tan 771 close to k = 0. From the definition in Eq. (|25[) . 
— tan rji /k 3 converges to the scattering length for p par- 
tial wave as k — ► 0. Alternatively, — tan 771 jk 2 is linear 
as k — > and the slope is the scattering length. 

A. Cesium 

Consider the scattering of two Cesium atoms in the 
3 S U state . The form of interaction potential is [Tt} 

v{R) = \br*^w-(% + % + ^)uR). 

(31) 

where B = 0.0016, a = 5.53 and = 1.072. These are 
given in ref [l2| and f c {R) is the cut-off function and has 
the expression 

f c (R) = G(R - R c ) + B(R C - R) exp-^/*- 1 ) 2 , (32) 

where 0(. . .) is the step function which is equal to 1 or 0, 
depending on whether its argument is greater or less than 
zero. Here R c is the cut off parameter. For comparison, 
the total interatomic potential V(R) and van der Waal's 
potential are shown in Fig[3J The ratio of the two are 
also shown as an inset plot. The figure shows that V(R) 
approaches the asymptotic form at around R = 24, where 
the ratio of the two is 0.8. Beyond R = 30, the two are 
indistinguishable. Using this potential, the phase shift 
and scattering lengths are calculated for the 132 Cs and 
134 Cs fermionic isotopes. As test calculations, we eval- 
uated the s wave phase shifts for the 133 Cs isotope. It 
is to be mentioned that, the value of 7 is 41240.1. How- 
ever, in ref. [l2| the value of 7 is defined as 41200. Using 
the later, the results of the test calculations are in good 
agreement with that of ref. [l2j ■ 

From Eq which provides a rough estimate, the 

value of a\ is 137a.u and —108 a.u. for 132 Cs and 132 Cs 
respectively. The 771 calculated from the variable phase 
method, for different R c for a range of k close to zero are 
shown in Fig[U Our calculations show that, the value of 
771 at higher values of k arc sensitive to 7, the dependence 
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tion is integrated till 771 converges to the order of 10~ 6 , 
which is consistent with the choice of tolerance. For Cs, 
this requires integration up to radial distances of ~ 500. 
The large radial distance is necessary as the phase shift 
is an asymptotic property. In addition, for the p wave 
scattering the centrifugal potential, which has R~ 2 de- 
pendence, vanishes at a slower rate compared to the van 
der Waal's potential at large radial distances. This slows 
the convergences of the phase shift. A calculation of the s 
wave phase shift confirms this. The s wave phase shift r]o 
converges to the order of 10~ 6 when the phase equation 
is integrated to a radial distance of w 100. 



FIG. 3: The black and red curves are the plots of interatomic 
and the van der Waal's potential respectively for Cesium. In 
the inset ratio of van der Waals term to the actual potential 
is plotted. 
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FIG. 5: Cube of scattering length a?(fe) or -tan 771 /fc 3 for l:i4 Cs 
is shown as a function of k. The black, green, blue and cyan 
correspond to R c = 23.115, 23.165, 23.190 and 23.215 respec- 
tively. 
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FIG. 4: The left and right panel of the figure are the phase 
shift ?7i(fc) of 132 Cs and 134 Cs isotopes respectively. Each 
panel show the plots of 771 (fc) for different cut off radius R c 
of the potential. The, black, red, green, blue and cyan cor- 
respond to R c = 23.115, 23.140, 23.165, 23.190 and 23.215 
respectively. 



is nonlinear as the equation of phase function is a non- 
linear differential equation. Consequently, small change 
in 7 could result in significant change of 771 . The depen- 
dence 771 on 7 for larger difference is evident from FigHJ 
which shows 771 for 132 Cs and 134 Cs. The difference in 
the interatomic potential of the two isotopes is the mass, 
which manifests as unequal values of 7. There are two 
distinct differences in the results for the two isotopes: for 
the R c = 23.215 case 771 converges to 179.08 and 182.21 
for 132 Cs and 134 Cs respectively; and for the R c = 23.115 
case 771 converges to 182.22 and 185.35 respectively. The 
phase function in the neighborhood of fc = is also sen- 
sitive to the accuracy of the integration. 

In the present calculations, the phase function equa- 



The values of — tan(?7 1 (fc))/fc 3 for 134 Cs within a range 
of fc for different R c are shown in Fig [51 As discussed 
in ref. [l2|, the R c = 23.165 is a physically reasonable 
choice. For this the scattering length calculated as the 
slope of — tan(77i)/fc 2 in the fc — > are 53 and —113 in 
atomic units for 132 Cs and 134 Cs respectively. This dif- 
fers from the analytic values mentioned earlier by 61% 
and 4.4% for 132 Cs and 134 Cs respectively. It shows that 
the analytic expression though crude provides a good es- 
timate for 134 Cs, which has higher 7 . 

The evaluation of 771 from the numerically calculated 
A/B poses some difficulty. This is to do with the choice 
of appropriate R*, the matching point, the dependence 
is shown in the Fig [6] The over all trend of the variation 
is rather complicated for the actual interatomic potential 
but less when only the van der Waal's potential is con- 
sidered. However, it is possible to calculate A/B from 
the value of a\ calculated earlier, which is obtained from 
the variable phase method. The value of A/B obtained 
from such a calculation are « —0.109 and « 1.042 for 
132 Cs and 134 Cs respectively. The value of R* around 
which A/B is close to this value occurs around the radial 
range of 23 — 24 atomic units. This is not surprising, 
as mentioned earlier, this is the radial range where the 
interatomic potential approaches the asymptotic form. 
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FIG. 6: \A/B\ for 134 Cs evaluated from the WKB solution is 
shown as a function of the matching point R* . The black and 
blue portions of the curve indicate the positive and negative 
values of A/B. The dashed red curve shows the A/B when 
only the van der Waal's penitential is considered and the inset 
plot shows the overall nature of the same. 



B. Lithium 

For the a 3 S+ state of 6 Li2, we use the interatomic po- 
tential suggested by Zemke and Stwalley [3. The long 
range part of the interaction potential, for 7 Li in partic- 
ular, are discussed in ref. [3]. For the R < 6.388674 and 
R > 18.0 regions, the analytic expressions recommended 
in ref. [3 are used. Then, for the 6.388674 < R < 18.0 
region, the interatomic potential is approximated as the 
cubic spline fitted to the potential values given by Zemke 
and Stwalley. A plot of the interatomic potential ob- 
tained is shown in Fig|7] As to be expected, the potential 
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K 



FIG. 7: Interatomic potential for aE„ of LI2. The solid 
(black) and dash (red) are the interatomic potential and the 
van der Waal's potential respectively. The plot in the inset is 
the ratio of the van der Waal's potential and the interatomic 
potential. 



approaches the asymptotic form, van der Waal's poten- 
tial, at smaller radial distance compared to Cs. This is 
evident from a visual comparison of the two potentials 
shown in Fig[7] and Fig[3] This implies that, the phase 
shift calculation of lithium can have better convergence 
properties. Consequently, the methods adopted for the 
Cs calculations are applicable to Li and provide results 
with higher accuracy. 

To optimize the calculations, we use the RKF method 
to solve Eq. (f2"8")) with relative tolerance varied as a func- 
tion of k. The relative tolerance is set to 10~ 5 for 
k > 0.01 and it is lowered as k is decreased. This is 
essential as the scattering length depends on the nature 
of the phase shift in the neighborhood of k = 0. Based on 
extensive test calculations, the optimal choice of the rel- 
ative tolerance is 10~ 6 for 0.01 < k < 0.005 and 10~ 7 for 
k < 0.005. With this choice it is possible to get reliable 
results for phase shift up to k = 0.0005. We find that 
phase shift stabilizes to 31.46 as k — ► which is shown 
in FigE 

33.0 



32.5 - 



* 32.0 - 



31.5 - 




31.0 

0.000 0.010 0.020 0.030 

k 

FIG. 8: The p-partial wave phase shift 771 (k) for the 3 Ej 
interaction potential of 6 Li2. The phase shift converges to 
31.46 as k -> 0. 

The sensitivity of the phase shift to the relative toler- 
ance at low k is not very prominent from the values of rji . 
However, it is not so with the — tan(?7i)/fc 3 , which is the 
scattering length in the k — > limit. The small value of k 
and the transcendental function tan enhances variations 
in 771. This is evident when — tan(?7i)/fc 3 as a function 
of k is compared at low k values for two different rela- 
tive tolerances. The plot of such a comparison is shown 
as the inset plot in FigO With improved tolerance, the 
lowest k value up to which the reliable phase shift can 
be calculated is decreased. There is however a limita- 
tion to decreasing the value of k, up to which the phase 
shift is calculated, by decreasing the value of relative tol- 
erance. Below a certain value of relative tolerance, the 
RKF method fails to perform the integration of the phase 
equation up to the cut off point. From the calculations, 
we find that — tan(r7 1 )//c 3 stabilizes to sa —92660 , this is 
evident from the plot in Fig|5]). The cube root of this is 
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FIG. 9: Variation of — tan(r;i)/fc 3 is plotted as a function of 
k. At low k values, the relative tolerance of the calculation is 
1CF 7 . The inset plot shows the sensitivity of — tan(?7i)/fc 3 to 
the relative tolerance. The red and black curves in the inset 
correspond to tolerance of 10 -6 and 10 -7 respectively. 



fitting of straight line over the data points. The slope of 
the least square fitted line is —92600, the cube root of 
which is —45 and can be considered as a reliable value of 
the p-wave scattering length. If we compare it with the 
analytic approximation of a± obtained by using Eq , (|24[) 
as the value of A/B, we find that percentage error in- 
curred by using analytic approximation is 43.74 percent. 
VI. CONCLUSIONS 



the a±. Another equivalent way to calculate the scatter- 
ing length is to evaluate the slope of — tan(f7i)/fc 2 in the 
k — ► limit. We calculate this by least square fitting the 
data points for — tan(^i)/fc 2 between 0.005 < k < 0.0005 
with an additional point corresponding to zero phase shift 
at k = 0. As is evident from Fig, (fit))) there is very good 



-o.oi 



FIG. 10: Plot shows least square fitted line (dashed red curve 
) on a{{k)k values in the domain 0.0005 < k < 0.005 (black 
crosses ). An additional point coinciding with origin is also 
considered while drawing the least square fitted curve. The 
slope of this line is ~ —92600. 



For p partial waves, the centrifugal term in the effective 
interatomic potential leads to a WKB solution, which is 
more complicated than the one without it as in s partial 
wave case. Consequently, the constant A/B of the par- 
tial wave solution evaluated has radial dependence and 
apriori. it is not possible to choose an appropriate R at 
which a reliable A/B can be evaluated. The approximate 
analytic expression to calculate A/B provide good results 
for 134 Cs but has large errors for 132 Cs and 6 Li. Based 
on our calculations, using the variable phase method, we 
estimate the p-w&ve scattering length of 6 Li is —45. 
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